Are the pentaquark sum rules reliable?* 
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We rewiew and scrutinize the existing mass determinations of the pentaquarks from the exponential Laplace Sum Rules 
(LSR). We do not find any sum rule window for extracting optimal and reliable results from the LSR, due to the unusual 
slow convergence of the OPE and to the exceptional important role of the QCD continuum into the spectral function in this 
channel. Instead, we use in this channel, for the first time, Finite Energy Sum Rules (FESR), which exhibit a nice stability 
in the QCD continuum threshold tc, at which one can extract, with a good accuracy, the mass of the lowest resonance. 
Including the D — 7, 9 condensate contributions in the OPE, we obtain Mb = (1513 ± 114) MeV, and the corresponding 
residue A|, ~ —(0.14 ^ 0.49) x 10~^ GeV^^, which favours the / = 0, J — 1/2, and negative parity S-wave interpretation of 
the ©(1540). However, our analysis indicates a degeneracy between the unmixed / = and I — 1 S-wave states. In the / = 0, 
J = 1/2, P-wave channel, we obtain, for the P-resonance, Mp ~ (1.99 ± 0.19) GeV and Ap ^ -(0.7 - 7.1) x 10"^ GeV", 
which we expect to be discovered experimentally. Our results also suggest that some intuitive choices of the continuum 
threshold used in the LSR literature are inconsistent with the FESR results. Finally, a study of the Q-K-N coupling using 
a vertex sum rule shows that, for the / = 0, S-wave channel, the leading OPE contributions only start to order in the 
chiral limit nis — 0, indicating that the O is very narrow. 



1. INTRODUCTION 

Recent experimental discovery of the 0(1540) as a nar- 
row KN bound state in 7-nucleon and 7-nucleus pro- 
cesses, e+e" and hadronic machines PP'' have stim- 
ulated renewed theoretical interests in hadron spec- 
troscopy |3i4j . In this paper, we shall critically re- 
analyze the mass determinations of the isoscalar / = 0, 
Q pentaquark mass from the exponential Laplace sum 
rules (LSR) 5 7 8 within the diquark scenario (3j and 
propose new analysis using Finite Energy Sum Rule 
(FESR). 

2. THE PENTAQUARK CURRENTS 

The basic ingredients in the resonance mass determi- 
nations from QCD spectral sum rules |9I10| as well 
as from lattice QCD calculations are the choice of 
the interpolating currents for describing the resonance 
states. Contrary to the ordinary mesons, where the 
form of the current comes from first principles, there 
are different choices of the pentaquark currents in the 
literature. The following analysis also postulates the 
existence of a strongly bound pentaquark resonance. 
We shall list below some possible operators describing 
the isoscalar / = and J — 1/2 channel^, which would 
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^However, this narrow state has not been confirmed by several 
experiments with high statistics and very good particle identifi- 
cation, using either e+e~ |3] and hadronic initial states. 
"^The isovector 1 = 1 current for S-wave resonance have been 
proposed by 1111 . We have also checked (see also 1121 ) that the 
tensor diquark current used in 1131 is an isovector 7 = 1 instead 



correspond to the experimental candidate 8(1540) Pp. 
Defining the pseudoscalar (ps) and scalar (s) diquark 
interpolating fields as: 

Qlli^) = [ulix)Cd,ix)] , 

= K(x)C754(x)] , (1) 

where a, 5, c are colour indices and C denotes the 
charge conjugation matrix, the lowest dimension cur- 
rent built by two diquarks and one anti-quark describ- 
ing the 8 as a / = 0, = 1/2+ S'-wave resonance is 
(see also 0)^: 

4 = e'^'^e'^^fe^f'^QZQ^^^Cs^ , (2) 
and the one with one diquark and three quarks is 0: 

ri^ = 7f ^""'^^fe i^eSeil^dc -{u^ d)} . (3) 

This later choice can be interesting if the instanton 
repulsive force arguments JS] against the existence of a 
pseudoscalar diquark bound state apply. Alternatively, 
a description of the 8(1.54) as a / = 0, = 1/2+ P- 
wave resonance has been proposed by [2 and used by 
|S] in the sum rule analysis: 

iD^Qlb)Qccih57f.Cs^ . (4) 

We have generalized this current by considering its 
mixing with the following one having the same dimen- 

of an isoscalar / = as confirmed by the authors in the revised 
version of their paper. In the following, we shall neglect the 
isospin breaking discussed in I14|, which would come from higher 
order diagrams in our analysis. 

^A negative parity state can be obtained by multiplying by 75 
the diquark operator. 



sion and quantum numbers: 

rit^ = ^'-'"'^'''^^'f'QZQlel.iD^^Csl) . (5) 

3. THE QCD SPECTRAL FUNCTIONS 

For the QCD spectral sum rules analysis, we shall work 
here with the two-point correlators: 



i / d'^x (0|r77^(:E)77^(0)|0), (6) 



built from the previous rj currents. It possesses the 
Lorentz decomposition: 



(7) 



• The QCD expression of the correlators associated to 
different choices of the currents is known in the litera- 
ture |5I7I8| to leading order of PT series and including 
the three first non-perturbative condensate {D < 5,6) 
contributions. We have checked the QCD expressions 
given there and agree with their results. However, at 
that approximation, we have added some missing con- 
tributions in j^. 

• We have not included in the OPE the contribution 
of the D = 2 tachyonic gluon mass induced by the re- 
summation of the PT series ^B] bearing in mind that 
this effect will be negligible, to the accuracy we are 
working, as illustrated in some examples jl7j . 

• We have included the D = 7,9 contributions into the 
QCD expression of the spectral function associated to 
the current in Eq. (0) , which we shall extensively study 
as a prototype example in this paper. In doing this cal- 
culation, we have worked in the chiral limit = 0, 
such that for consistencies, we shall use the SU{3) 
symmetric value (ss) — (dd) for these contributions. 
In this particular example, we have checked that the 
contribution of the four-quark condensate vanishes to 
leading order. In our preliminary results, we also found 
that its radiative correction though not identically zero 
gives a negligible contribution. We have also neglected 
the contributions of the three-gluon condensate of the 
type g{GGG) (ss) assuming that the theorem in 18 for 
the light quark bilinear operators continues to hold for 
the diquark correlators ^, which factorize during the 
evaluation of the QCD expression. 

• We have evaluated the new contributions associated 
to the current rjnew in Eq. jSJ, where we found that, 
to leading order in and in the chiral limit niq — > 0, 
the contribution to the correlator vanishes. This result 
justifies a posteriori the unique choice of operator for 
the P-wave state used in |3I8| . 

4. THE LAPLACE SUM RULES (LSR) 

We shall be concerned with the Laplace transform sum 
rules: 

r°° 1 
£^/^{t) EE / dte-'^ -lmA"/B"{t), 

' Jt< 



(8) 



®We plan to check explicitly this result in a future publication. 



where i< is the hadronic threshold, and H denotes the 
corresponding hadron. The latter sum rule, or its slight 
modification, is useful, as it is equal to the resonance 
mass squared, in the usual duality ansatz parametriza- 
tion of the spectral function: 

- W^/i?^(i) ~ {Xl/XlMH)5{t - Mfj) + 

TT 

"QCD continumii" e(i - Q, (9) 

where the "QCD continuum comes from the disconti- 
nuity of the QCD diagrams, which is expected to give a 
good smearing of the different radial excitations . Xh 
is the residue of the hadron H; tc is the QCD con- 
tinuum threshold, which is, like the sum rule variable 
r, an (a priori) arbitrary parameter. In this paper, we 
shall look for the r- and ^c-stability criteria for extract- 
ing the optimal results. For illustrating our analysis, 
we give below the checked and completed LSR of the 
S'-wave current in Eq. ^ including the new D — 7,9 
high-dimension condensates in B: 
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pn being the notation in ^HIj while: (ss), (a^G^) are 
respectively the dimension D = 3 quark and D — A 
gluon condensates; g{saGs) = g{sa^'' {Xa/'2.)G'^jjS) = 
Mq{ss) is the D — 5 mixed condensate. Throughout 
this paper we shall use the values of the QCD param- 
eters given in Table 1. 

Table 1 

QCD input parameters used in the analysis. 



Parameters 



References 



fnj{2 GeV) = (111 ± 22) MeV .10J9i20l2ll 
(Jd)i/3_(2 GeV)=-(243± 14) MeV |10I19I22| 
(ss)/(dd) ==0.8±0.1 |l()ll9l23| 
{a,G^) = (0.07 ± 0.01) GeV^ 
Mg = (0.8±0.1) GeV^ fTTl^ 

We study the LSR in Eqs. (HUJ and ((TTJ- We find that 
all LSR corresponding to different currents present the 
common features shown in Fig. ^ 

— The _B-component increases rapidly with r. 
Then, it is useless at that approximation of the 

^A, the mass prediction decreases smoothly 
when T increases. The OPE converges for r < 0.9 
GeV~^ (LHS of the vertical dashed line). 
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Fi gurel. T-behaviour of A/@ for given values of tc- OntheLHS 
of the vertical dashed line, the OPE converges. The vertical line 
with arrow on each curve shows that the continuum contribution 
dominates over the resonance in the LHS region. A (resp B) 
corresponds to the invariant defined in Eq. 

— The QCD continuum contribution dominates 
over the resonance one in all ranges of r where 
the OPE converges (LHS of the vertical line with 
arrow on each curve). Indeed, for r < 0.9 GeV^^, 
the QCD continuum contribution eats more than 
84% of the OPE one. 

• Therefore, it is impossible to find a sum rule win- 
dow region where both the resonance dominates over 
the QCD continuum contribution, and where the OPE 
converges. Intuitively, this feature is expected as the 
current describing the pentaquark is of higher dimen- 
sions, and therefore is more affected by the continuum 
contribution than the well-known sum rule for ordinary 
qq mesons. The absence of the sum rule window is re- 
flected by the increase of the mass predictions with the 
QCD continuum threshold tc- In existing literature, 
the tc-values have been fixed ad hoc and intuitively. 

• During the evaluation of the different QCD diagrams, 
we do not find (to leading order in as) any factoriza- 
tion of the {su)-{udd) diagram corresponding to a re- 
ducible K-N continuum diagram, which has nothing 
to do with the diquark picture. Then, our direct ob- 
servation does not support the criticisms raised in [221 
and refuted in [23 on a possible double counting due 
to the non-subtraction of the reducible diagram in the 
existing sum rules analysis of the O. 

• We conclude from the previous prototype example 
that the LSR using the simple duality ansatz: reso- 
nance-|-QCD continuum criterion is not appropriate for 
determining the pentaquark masses due to the absence 
of the usual sum rule window. Due to the huge con- 
tinuum contribution (« 85%) at relatively large r « 1 
GeV~^, the LSR cannot strictly indicates the existence 
of the resonance into the spectral function. 

• We have checked (though not explicitly shown in the 
paper) that the conclusions reached in the paper also 
apply to the sum rules used in the literature: |2] (cur- 
rent in Eq. Q) and susbsequent uses in j2(jl27| for 
the / = 0, S'-wave state; the sum rules used in for 
the / = 1, S'-wave state; in |S] (current in Eq. 0)) 
for the / = 0, 1 P-wave state; in for the 1 = 1 
tensor current and the sum rules used in |28l29j for 
studying the J = 3/2 states. Indeed, in most LSR, the 
OPE does not converge at the scale where the results 



are extracted, while the QCD continuum threshold has 
been taken arbitrarily or intuitively. 
• The above results raise some doubts on the valid- 
ity of the results obtained so far in the existing litera- 
tures. Indeed, if one insists on using the LSR for pre- 
dicting the parameters and some other pentaquark 
states, it is mandatory to introduce a more involved 
parametrization of the continuum spectral function. 

5. FINITE ENERGY SUM RULES (FESR) 

In contrast to the LSR, Finite Energy Sum Rules 
(FESR) pOIHlllOj have the advantage to project out 
a set of constraints for operators of given dimensions 
(local duality) . They also correlate the resonance mass 
and residue to the QCD continuum threshold so 
avoiding inconsistencies of the values of these param- 
eters. Also contrary to the LSR, the resonance and 
QCD continuum contributions are separated from the 
very beginning. The FESR read: 

Mi^iA/B) = r dteimA" /B"\exp 



ta 



t< 



dt e linA"lB"\QCD 



(13) 



From the expressions of the spectral function given pre- 
viously, one can easily derive the FESR constraints. 
Doing the FESR analysis for the A{q^) invariant, one 
can notice that, at the approximation where the OPE 
is known (£* < 6), one does not have a stability in tc for 
different moments Jvi^ and for different choices of the 
currents (see Fig. Therefore, we will not consider 
this invariant in the paper. 

The / = 0, S-wave channel 

We illustrate the analysis by the current in Ref. [S] 
(the other choice [2| in Eq. JSJ) has approximately the 
same dynamics as one can inspect from the QCD ex- 
pressions) . Including the D = 7 and 9 condensate con- 
tributions, the two first lowest dimension constraints 
from the B{q^) invariant read: 
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from which one can deduce the mass squared: 
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Figure 2. tc-behaviour of Mq for different truncations of the 
OPE. 

The behaviour of Mq is shown in Fig. 13 for different 
truncations of the OPE. One can notice a stabihty at 
tc — 2.29 GeV^ , where the OPE starts to converge 
after only the inclusion of the D — 7 + 9 condensates, 
while D — 7 alone destroys the stability reached for 
D < 5. One can notice the important contribution of 
the lowest quark and mixed quark-gluon condensates 
in the OPE, which play a crucial role in this mass de- 
termination. To that order, we obtain: 

Me ~ (1513±20±10±40±30±30±95) GeV 

~ (1513 ± 114) GeV , (17) 

where the errors come respectively from 
nis, {qq), {dsG"^), Mq, the estimate of the higher di- 
mension condensates and the violation of the vacuum 
saturation assumption of the D — 7,9 condensates by 
a factor (2 ± 1) like in the p-meson |32I3()| and some 
other channels ^Ol- One can notice that: 

• The existence of the <c-stability point makes the 
superiority of FESR compared to the LSR in this 
channel. For the LSR, Mq increases with tc- Here, 
the localisation of the stability point induces here a 
negligible error. 

• The FESR order parameter in the OPE, tc ~ 2.3 
GeV^ is much larger than for the LSR (r^^ < 1 GeV^), 
implying a much better convergence of the OPE for 
the FESR, and then a much more reliable result than 
the LSR. 

• Working with ratio of higher moments 
A^f .leads to almost the same value of Mq. 
The slight variation is much smaller than the error in 
Eq. (dZj). 

• Truncating the OPE at D = 5 like done in the avail- 
able literature would give a slightly lower value of Mq 
at the stability point tc ~ 3.2 GcV^ (see Fig. but 
compatible with the one in Eq. (|17|l . 

• Contrary to Mq, the value and the sign of Aq are 
very sensitive to the truncation of the OPE due to 
the alternate signs of the condensate contributions in 
the analysis. The stability in tc is obtained after the 
inclusion of the D = 5 condensate contributions as 
shown in Fig. |21 To our approximation D < 9, the 
most conservative result is: 



Figure 3. tc-behaviour of A 
condensates in the OPE. 



10'' in GeV'^ including D < 5 



where the range comes from the shift of the fc-stability 
point from D = 5to D = 9 approximation. This result, 
though inaccurate, suggests that the parity of the O is 
negative, as indicated by the lattice results given in 
|33|^. Improving the accuracy of our result requires 
more high-dimension condensate terms in the OPE. 



The I ~ 1, S-wacve channel 

We have also applied FESR in the / 

nel with the current |11| : 

1 --r^. 



1 S'-wave chan- 



V2 



where t is an arbitrary mixing parameter. To the 
D < 5 approximation, the analysis gives almost the 
value of Ajfe hi Fig. |21at the same approximation. This 
result can be interpreted as a consequence of the good 
realization of the SU{2)p symmetry for the u and d 
quarks. Then, we expect that the unmixed / = 1 part- 
ners of the unmixed / = state will be around the 1.5 
region if any. 

The / = 0, _P-wave channel 

We do a similar analysis for the P-wave current given 
in Eqs. Q and (jSJ, where as we have mentioned in 
section 3, the contribution from Eq. (|SJ) vanishes to 
leading order in as- The corresponding FESR up to 
D — 5 condensates are given below ^: 
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To this order, the moments present a similar tc- 
behaviour as above (see Figs. 0]and[Sl). Both for 
the mass and residue, the stability point is: 

tc ~ 5.5 GeV^ 



(22) 
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^At the approximation D < 5 the LSR does not converge such 
that analogous results obtained in I5l26l27i should be taken with 
a great care. 

**The inclusion of higher dimension condensates is in progress. 
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Figure 4. tc-behaviour of Mp including the D < 5 condensates 
in the OPE. 

The corresponding resonance mass and residue are: 
Mp ~ 1.99 ± 0.19 GeV, 

Xl w -(0.7 - 7.1) X 10"^ GeV" . (23) 

The errors come mainly from the estimate of the un- 
known D = 7, 9 condensates contributions inspired 
from the S-wave channel. One can notice that: 
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Figure 5. tc-behaviour of A|, X 10^ GeV^^ including the D < 5 
condensates in the OPE. 

• The value of the QCD continuum threshold at which 
the FESR stabihzes is much higher than the intuitive 
choice used in the LSR IH^ needed to reproduce the ex- 
perimental mass of the O. 

• The mass value obtained for the P-wave resonance is 
(450 ± 190) MeV higher than the 9(1540) mass, which 
suggests that there is a P-wave state different from the 
0(1540) in the region around 2 GeV, which we expect 
to be discovered experimentally. 

• The value and sign of the residue suggest that this 
P-wave state has a negative parity like the 0(1540). 

6. THE e-K-N COUPLING 

For studying this coupling, we start from the three- 
point function: 

V{p,q) ^i^ j d*x ct^y e*(P^+9!') {{)\ri{Q)N{x)K{y)\Q) (24) 

where K{y) = {rus + mu)s{ij5)u is the kaon current, 
while N(x) =: u{C^z)du : +b : u[Cd)"f5U : is the nu- 



cleon interpolating field |34j {b being an arbitrary mix- 
ing parameter), rj is the O current defined in previous 
section. For definiteness, we work with the 5'-wave 
current given by A QCD evaluation of the vertex 
in the chiral limit nig = shows that the leading and 
Us orders perturbative and non-perturbative diagrams 
give zero contributions ^ . The result then suggests 
that the Q-K-N coupling is of the order supporting 
the experimental observation ^ that the 0(1540) is a 
narrow state. The narrowness of a pentaquark state 
has been already advocated in the past, from duality 
arguments, where its decay into BBB baryon states is 
dynamically favoured, implying that light pentaquark 
states below this threshold are naturally narrow |14 |. 
A narrow S-wave pentaquark state has been also ob- 
tained in j2S| and ^3] using simple chiral symmetry 
arguments. In the narrowness of the O is due to a 
destructive interference between two almost degenerate 
states, while in |37], it is due to the flavour structure 
of the O, which, after the meson formation, the resid- 
ual three-quark piece has a little overlap with the octet 
baryon wave- function. 

7. SUMMARY AND CONCLUSIONS 

• We have re-analyzed the existing LSR results in the 
literature. We found that due to the slow convergence 
of the OPE and to the relative importance of the QCD 
continuum contribution into the spectral function, the 
minimal duality ansatz "one resonance + QCD con- 
tinuum" is not sufRcient for finding a sum rule window 
where the results are optimal. These features penal- 
ize all existing sum rule results in the literature, which 
then become unreliable despite the fact that the mass 
predictions reproduce quite well the data. However, 
this apparent good prediction is due to the fact that the 
intuitive or arbitrary choice of the continuum threshold 
tc- In fact, in the LSR analysis, the mass prediction 
increases with though it is a smooth function of r 
as can be seen in Fig. ^ 

• On the contrary, FESR has the advantage to present 
a good tc-stability and converges faster than the LSR, 
because the optimal results are obtained at higher scale 
tc~(2^ 3) GeV^ than the one of LSR t'^ < 1 GeV^. 

• Truncating the OPE at the D — 9 condensates, at 
which the OPE starts to converge, we obtain the re- 
sult in Eq. H17|) . for the S'-wave state, which one can 
compare with the experimental candidate 0(1540). 

• By truncating the OPE at D = 5, we also find from 
FESR a good degeneracy between the unmixed / = 
and I — I S'-wave states. 

• Similarly, we obtain, from FESR, the mass of the 
P-wave state of about 2 GeV in Eq. H23|l to order 
£> = 5 of the OPE, but including the estimated effects 
of £> = 7, 9 condensates. This mass is (450 ± 190) MeV 
higher than the 0(1540). 

• Finally, an analysis of the Q-K-N coupling using ver- 
tex sum rules supports results in the literature that the 

^Our result is stronger than the one in Ref. 13 which claims a 
non-zero as contribution. 



0(1540) is a narrow state. 

• Our results seem to favour the case (b) discussed 
in where the Q resonance is induced in KN scat- 
tering by coupUng to a confined channel. A complete 
program using FESR in different pentaquark channels 
is in progress and will be published elsewhere. 
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